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Abstract 

We consider the vacuum structure of the finite M = 2 theory with Nf = 2N fundamen- 
tal hypermultiplets broken to M = 1 by a superpotential for the adjoint chiral multiplet. 
We do this in two ways: firstly, by compactification to three dimensions, in which case 
the effective superpotential is the Hamiltonian of an integrable spin chain. In the second 
approach, we consider the Dijkgraaf-Vafa holomorphic matrix model. We prove that the 
two approaches agree as long as the couplings of the two theories are related in a partic- 
ular way involving an infinite series of instanton terms. The case of gauge group SU(2) 
with Nf = 4 is considered in greater detail. 



*e-mail: t . hollowood@swansea .ac.uk 

'e-mail: k.ohta@swansea.ac.uk, k-ohta@riken.jp 



1 Introduction 



In this work, we examine the breaking of the M = 2 SU(N) theory with Nf = 2N massive 
fundamental hypermultiplets to M = 1 by adding a general tree-level superpotential for 
the adjoint field of the form 



Since we describe the situation with general hypermultiplet masses, the same problem 
with Nf < 2N can be described by taking limits of large masses where the associated 
hypermultiplets decouple. 

There are several known ways to tackle the problem: factorizing the Seiberg-Witten 
curve [T]; using MQCD or deformation of the Seiberg-Witten curve jSJEj; and using matrix 
models and the glueball superpotential (for a review see jl] and references therein). As 
well as describing the matrix model technique in some detail, we will add another approach 
to the list based on integrable systems. 

The integrable system approach begins with the observation that the Seiberg-Witten 
curve of a class of M = 2 theories is the spectral curve of an integrable system. If one 
compactifies the theory on a circle to three dimensions and breaks to M = 1 as in (jl.lj) . 
then one can write an effective superpotential for a set of scalar fields that includes the 
scalar from four dimensions as well as the dual photons and Wilson lines of the gauge 
field. The Coulomb branch in three dimensions is identified with the (complexified) phase 
space of the integrable system and the effective superpotential is one of the Hamiltonians 
of the integrable system. Hence, vacua correspond to equilibria of the integrable system 
(2]. The cases which have been considered in the literature are: 

(i) The pure M = 2 theory in which case the integrable system is the affine Toda 
system 00 E|. 

(ii) The M = 2* theory (massive adjoint hypermultiplet) in which case the integrable 
system is the elliptic Calogero-Moser system [HUSJ- 

(iii) The M = 2 finite quiver theories in which case the integrable system is the spin 
elliptic Calogero-Moser system [TO] . 
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(iv) The compactified five-dimensional M = 2* theory in which case the integrable 
system is the elliptic Ruijsenaars- Schneider system [TTj . 

(v) The Leigh-Strassler deformed M = 4 theory in which case the integrable system is 
also the elliptic Ruijsenaars-Schneider system ^21 G2] (i n a sense this example is somewhat 
different because there is no M = 2 parent theory). 

In fact the integrable system approach can be adopted in an abstract sense by con- 
sidering the problem at the level of action/angle variables [S] even when no concrete 
realization of the integrable is known or can easily be written down. Now we turn to the 
theory with fundamental hypermultiplets. The M = 2 theory with vanishing beta func- 
tion, Nf = 2N, is known to be associated to a spin-chain integrable system jTU UHl EH Ej 
and this suggests that we can associate the vacua after breaking to Af = 1 to equilibria 
of this integrable system. In Section 13.21 we consider the relation of the Nf = 2N theory 
with the spin chain integrable system in more detail and use it to deduce the vacuum 
structure of the SU (2) theory with four massive hypermultiplets. 

We note that the formalism can in principle be extended beyond SU(2) but it becomes 
harder to find all the equibilibria. In Section we use the alternative matrix model 
technique to find the glueball superpotentials of these theories, an approach that readily 
extends to Nf < 2N. We note that the coupling of the matrix model approach is not 
equal to the coupling of the Seiberg-Witten curve. The two are related by a series of 
instanton terms: 

oo 

ri = r 2 + J2c k e mkT2 . (1.2) 

fc=0 

The actual relation will be determined in Section 5. In fact it is known that the coupling 
of the Seiberg-Witten curve is itself not equal to the bare coupling of the theory j!8j . 
rather they are also related by an expansion of the form (|1.2jl . 

2 Semi-Classical Analysis of SQCD Vacua 

The tree level superpotential of A/"=l SU(N) supersymmetric QCD with Nf flavors is 
obtained from the J\f=2 superpotential by the addition of a supersymmetry breaking 
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potential for the adjoint scalar field: 

W tree ($,Q,Q) = TrQ$Q-TrMQQ + Tr W{$) , (2.1) 

where Q and Q are NxNf and Nf x N matrices, respectively, and M = diag(mi, . . . , mjvy) 
is the mass matrix. In what follows it will be convenient to formulate the SU(N) theory 
in terms of the U(N) theory with an explicit Lagrange multiplier constraint in order to 
enforce the tracelessness of $. To this end we take the potential for the adjoint to have 
the form 

n+l 

W(x)=^ii i x i -£x, (2.2) 
i=i 

where £ is the Lagrange multiplier. 

To investigate the vacuum structure one imposes D- and F-flatness modulo gauge 
transformations. The latter conditions read 

QQ + W\$) = , $Q = QM , Q$ = MQ . (2.3) 

As usual solving the D- and F-flatness conditions modulo gauge transformations is equiv- 
alent to only imposing F-flatness modulo complex gauge transformations. The latter can 
be used to diagonalize $: 

$ = diag (0i,0 2 ,..., 0^) . (2.4) 

This leaves the abelian U(1) N subgroup as well as the elements of the Weyl group which 
act as permutations of the elements ()2.4j) unfixed. For generic masses the solutions are as 
follows. First of all, for r = 0, . . . , min(iV, Nf), 

$ = diag (m 7l , . . . ,m 7r ,0 r+ i, . . . ,0at) , (2.5) 

where the Ij are distinct elements of the set {1, . . . , Nf}. The <pj, j = r+1, . . . ,N satisfy 

W'^> j )=0, j=r + l,...,N. (2.6) 

= dj6 n . , Qij = dj6 n . , j = 1, . . . , r , (2.7) 
and otherwise, and where 

djd, + W\m h ) = , j = 1, 2, . . . , r . (2.8) 

The fixed gauge symmetry and remaining complex gauge symmetry can be used to set 
di — 1 and so the solution is at least generically discrete. 
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The solution breaks the gauge symmetry to Yli=i U(Ni) with J2i=i N% = N — r. 1 Here, 
Ni is the number 4> a which lie at the i th solution of the polynomial W'(x) = 0. If W* mc 
is the effective tree-level superpotential in the vacuum then variable £ is then determined 
by setting dW*Jd£ = 0. 

Now consider the situation with the simplest potential W(x) = fix 2 /2 — £x. In this 
case, 

& = £A*, j = r + l,...,N. (2.9) 
The effective superpotential is 

W.~ = E(H-^)-^ 2 - (2-10) 

Hence, 

j=i 

and so in the SU (N) theory we have 

K«=fE< + ^)(E^) 2 - < 2 - 12 > 

j=l V ^ 3=1 

In this vacuum SU(N) is broken to SU(N — r) and for r < A/" the vacuum is characterized 
by the meson VEV 

M = QQ = diag (di, . . . , d r , 0, . . . , 0) . (2.13) 

In the low energy limit, the theory is an M = 1 theory with SU(N — r) gauge group 
which conventional reasoning says will confine and there will be iV — r vacua. Therefore 
we expect the total number of these vacua is 

min(JV,JV»-l 

E (N-r) Nf C r , (2.14) 

r=0 

where N f C r represents the number of subsets {ii, . . . , I r } C {1, . . . , Nf}. 

The case with r = N, which means Nf > N, is somewhat special but can only be 
attained if the iV-dimensional subset of the masses satisfy the constraint YliLi m ij = 0- 
1 Note that the overall U(l) is frozen out and non-dynamical. 
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In this case the vacuum is characterizes by a non-zero baryon VEV: 

B h - lN = ^ ■ j -Q hIl ■ ■ ■ Q jNlN = d x d 2 ■■■d N , 

1 . . . . ( 2 - 15 ) 

B h~iN = J^ e h-jNQiih ' ' ' Qinjn = d\d 2 •••dN . 

In these "Baryonic" vacua the gauge group is completely broken. 



2.1 Example: SU{2) with N f = 4 

In this section, we consider the finite theory with gauge group SU (2) and Nf = 4 with 
a quadratic superpotential W(x) = fix 2 /2, i.e. in our approach we take the U(2) theory 
with W(x) = [ix 2 /2 - £x. 

We start with the vacua with r = 0. In this case there is no Higgsing and 

$ = Q = Q = 0. (2.16) 

At low energy, the theory is described by a pure Af = 1 theory with 577(2) gauge group. 
This will confine and lead to 2 vacua. The superpotential vanishes in the classical ap- 
proximation but we expect it to receive corrections from fractional instantons with the 
characteristic behaviour e mT . 

There are four vacua with r = 1 with, for I = 1, . . . , 4, 

$ = diag (mj, tin) , i = -firm (2.17) 

For each choice I, Qu and Qn are non-vanishing. The SU(2) gauge group is completely 
broken by the Higgs mechanism and the classical value of the superpotential is 

W^ = iMmj. (2.18) 



The baryonic vacua can only occur when the masses are non-generic, in fact one needs 
m; + mj = 0, for some I and J. In this, one can check that there is a moduli space of 
vacua. 
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2.2 A brane configuration 



The vacuum structure that we have established on the basis of a classical analysis can 
be readily interpreted in terms of the brane configurations. The SU(N) gauge theory 
with Nf flavors that we are considering is realized by a configuration of N finite D4- 
branes stretching between two NS5-branes and Nf semi-infinite D4-branes emerging from 
one of the NS5-branes in Type IIA string theory ^Hj (see the review |2D] and references 
therein). The world- volumes of D4-branes and NS5-brane are along the 01236 and 012345 
directions, respectively, and the gauge theory arises as the effective theory on the common 
directions 0123. 

For a theory with J\f—2 supersymmetry, the two NS5-branes are parallel and the D4- 
branes corresponding to gauge and flavor symmetry are freely moving on 45-plane. The 
x 4 + ix 5 positions of finite D4-branes correspond to the (classical) VEV of the adjoint 
scalar field $ and ones of semi-infinite D4-branes correspond to the masses of the matters 
in fundamental representation. 

The supersymmetry breaking from Af=2 to A/"=l is achieved by a modification of the 
relative "shapes" of NS5-branes. For example, turning on the simple mass perturbation 
W(x) = \xx 2 /2 — £x, 2 is achieved by rotating one of the NS5-branes into the 78 direction 
[2*T] . Denoting w = x 7 + ix s , the the rotated NS5-brane is described by the complex line 
w — \ix — £. The other NS5-brane remains along w — 0. Generally, deforming the M—2 
theory by turning on the k + 1-th order potential as in (j2.2j) is achieved by demanding 
that the shape of one of the NS5-branes is described by w = W'(x). 

The relation between the classification of vacua as in Section |21 and the brane picture 
is rather direct. A vacuum with given value of r, corresponds to a situation in which 
r of the semi-infinite D4-branes match up in x-space at x — mi i i = 1, . . . ,r, with r 
of the finite D4-branes to make semi-infinite D4-branes which now end on the left-hand 
NS5-brane. The remaining N — r finite D4-branes have to lie at zeros of W'{x). 

For example, for gauge group SU{2) with Nf = 4, and the quadratic deformation 
W(x) = [ix 2 j 2 — £x, there are generically six vacua. There are 2 confining vacua where 
the two finite D4-branes lie at x = w = giving rise to an unbroken SU(2) in the IR 

2 Once again we introduce the Lagrange multiplier £ to freeze out the U(l) degree-of- freedom. The 
variable £ is fixed by demanding the centre-of-mass of the finite D4-branes is zero in the x plane. 



6 




(a) (b) (c) 



Figure 1: The classical brane configuration for (N,Nf) = (2,4). In general, the configu- 
ration has two finite D4-branes between the rotated NS5'-brane and NS5-brane, and four 
semi-infinite D4-branes from the right NS5-brane. The bullets on the line denote the end 
point of D4-branes on the NS5-brane. (a) Confining vacua. The number of vacua is 2. (b) 
Higgs vacua. There are 4 different choices in total, (c) "Baryonic" vacua for non-generic 
masses. The finite D4-branes disappear and brane configuration splits. 

and the usual count of 2 vacua. There are 4 Higgs vacua where one of the semi-infinite 
D4-branes links up with one of finite D4-branes at x — mi, I G {1,2, 3, 4}. The remaining 
finite D4-brane sits at x — mi. The confining and Higgs vacua are illustrated in Fig. 1 
(a) and (b). Non-generically there are other possibilities. For example, if mi + mj = 
then both finite D4-branes can link up with 2 of the semi-infinite D4-branes. In this case, 
as illustrated in Fig. 1 (c), the brane configuration splits into two separate pieces. 

3 Integrable Systems and Vacuum Structure 

3.1 The general story 

In this section we set out a general picture of the relation between TV = 2 theory broken 
to M = 1 and integrable systems. In this story it is very useful to have in mind the brane 
configurations of Type IIA string theory and brane rotation/deformation. 

It is well-established that J\f = 2 theories are related to integrable systems ^7j. The 
method that we describe in this section exploits this fact. The Seiberg-Witten curve 
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of the theory is identified with the spectral curve of the integrable system. The moduli 
space, or Coulomb branch, of the theory, is therefore identified with the space of conserved 
quantities in the integrable system. The Coulomb branch has a natural set of coordinates 
which are defined by integrating the Seiberg-Witten holomorphic 1-form A around a set 
of canonical A-cycles on the curve S: 

a l = I A . (3.1) 

Here, i = l,...,g, the genus of S. The cycles are chosen so that the ctj determine the 
masses of electrically changed BPS states int he theory. Magnetically charge BPS states 
involve the dual quantities 

aP = I A , (3.2) 



where the cycles Bi are chosen to have intersections Ai R Bj = Sij. Remarkably, the a; 
are precisely the action variables of the integrable system. The conjugate angle variables 
do not appear directly in the theory, however, they are identified with the coordinates ipi 
of a point in the Jacobian of E. The Jacobian torus is defined as follows (for a reference 
on Riemann surfaces see [22] )• Let LUj be the associated set of g holomorphic 1-forms 
(abelian differentials of the 1st kind) normalized so that § A u k = Sjk- The period matrix 
of S int is the g x g matrix with elements 



T jk = (h LU k . (3.3) 

Jb, 

The Jacobian torus consists of points ipj £ C 9 with the identifications 

ipj ~ ipj + rij + r jk m k , rij, m k £ Z . (3.4) 

As defined the Oj and ipi are canonically conjugate: 

{ai,ipj} = 5ij . (3.5) 

In terms of the action angle variables, the dynamics with respect to any choice of Hamil- 
tonian function H^oA is particularly simple: 

hi = , ipi = ^ . (3.6) 

We can be more explicit by using the fact that the spectral curve of the integrable 
system is described by a relation of the form F(t,x) = which is nothing but the usual 
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description of the Seiberg-Witten curve of the field theory coming from the Type IIA 
brane configurations described in Section 17121 In particular, the Seiberg-Witten 1-form is 



A = x - . (3.7) 

An important role is played by the points Pj, i = 1, . . . , k + 1, where x = oo at which E 
is non-compact. In the brane configuration these are the points where one of the k + 1 
NS5-branes goes to infinity. In order to see their significance in the integrable system, let 
us consider the flows generated by some Hamiltonian at the level of the angles variables 
(13. 6 j) . It turns out that each flow can be written in terms of a meromorphic 1-form Q on 
£ (actually an abelian differential of the 3rd kind) which is normalized by the condition 3 

tt = . (3.8) 
The angular velocities are then given by integrals around the dual cycles: 

r oh 

The 1-form Q has poles at the points Pj and the asymptotic expansions at these points 
are determined by the choice of Hamiltonian H: 

lim ft(p) = d(Ui(x(p)) + 0(l/x(p))) , (3.10) 

P-*Pi 

where the Ui(x) are polynomials which we will relate to the tree-level superpotentials 
which break M — 2 — ► M — 1. In fact in the brane picture w = Ui(x) describes the 
bending of the i th NS5-brane in the 89 direction asymptotically at infinity. The M = 1 
deformation of the i th gauge group factor is associated to a potential 

W' l {x) = U l+l {x)-U l {x) . (3.11) 

Note that a shift in the Ui(x) — > Ui(x) + f(x) is physically irrelevant. In the integrable 
system, the polynomials Ui(x) determine specify Q uniquely, up to exact terms which 
correspond to the freedom to shift each Ui(x) by the same function. In this sense, only 
the differences of the Ui(x) matter — as (|3.11j) makes clear. 

We now relate Q, modulo exact terms, to H. The first move is to employ the Riemann's 
bilinear relation to the abelian differentials 4 

d\ 

it, = — 3.12 

da; 



3 This can always be arranged by adding a suitable linear combination of the u>i 
4 Note that it will be important below that the a,i derivative is taken at fixed x. 
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and fi: 
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r~lb (<f n f f / ^) = J^Resp^Uii^j) , (3.13) 

i = \ \" Ai J Bi J Bi J Ai J i=l 

Now we use the normalization condition ()3.8j) and the fact that (J3.1)) and (J3.12)) imply 
that 

r tt, = Sa , (3.14) 



I 

J A, 



to deduce 

fc+l o fc+i 



# ft = V'Resp i (C/ i (a;)7r i ) = — Res Pi (£/i(x) logtdx) . 
i=i dQi i=i 



(3.15) 



In the last line we used the fact that the a, derivative is at constant x and so could be 
pulled outside the integral. Hence, by comparing this with (|3.fi|) we deduce that 

k 

H = -^Res Pi (Ui{x)\ogtdx) (3.16) 



At this point we can investigate the equilibrium configurations of the integrable system 
which ought, by our general reasoning, to be associated to vacua of the QFT. For an 
equilibrium, 

^ = / O = . (3.17) 

Odi JBi 

Given the normalization condition this implies that for an equilibrium Q = dw, for some 
meromorphic function w with singularities at Pi with asymptotic behaviour from (|3.10p . 
w — > Ui(x) H — • . Hence the existence of a vacuum of the gauge theory can be phrased as 
the existence of a certain meromorphic function on the curve of the theory. The relation to 
the brane configuration should now be obvious. The points Pi, as we have said, represent 
the asymptotic positions of the NS5-branes. The meromorphic function w describes the 
deformation of the Seiberg-Witten curve into the 89-plane associated with the breaking 
of AT = 2 to Af = 1. The superpotential ()3.1(i|) at the critical point can be written as an 
integral over the Riemann surface 

H* = — — / \ogtdwAdx, (3.18) 

from which one can relate to Witten's MQCD superpotential [23J (see l3j for a discussion 
of the relation). 
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For our particular example, there are two NS5-branes and hence two points at infinity 
Pi = oo + and P 2 = oo~. We can choose U 2 = and U\{x) = W'(x). By integrating-by- 
parts in (jH.lfij) . we have a Hamiltonian and hence an effective superpotential 



Notice that we have developed the theory of the integrable system entirely in the 
language of action/angle variables. However, for certain examples we can write down 
the Hamiltonians in terms of other kinds of variables; for instance, the positions and 
momenta of particles, or spin variables and mixtures thereof. We do this for the theory 
with Nf = 2N fundamental hypermultiplets in the next section. 

3.2 Spin Chains 

In this section, we describe how to extend the integrable system approach to describe 
the M = 2 theories with fundamental matter. The integrable system that is associated 
to the finite theory, with Nf = 2N, has been identified in EH EH El as a certain 
twisted 577(2) spin chain. This system describes the interactions of a set iV spins 
i — I, . . . , N, which are 3- vectors 




(3.19) 



° ~ \°1 i °2 i J 3 ) ■ 



(3.20) 



The overall length of each vector is non-dynamical: 




The parameters fc« will be associated to masses in due course. 



The integrable system has a 2 x 2 "Lax matrix" defined at each "site" : 




(3.22) 



The "transfer matrix" is then defined to be 



T(x) = Li(x) ■ ■ ■ L N (x)V 



(3.23) 
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and we have allowed for a "twisting" described by the constant matrix V. Unlike 
Uni E] , we will make the choice 




V = - , . J > ( 3 - 24 ) 



where h is a parameter that we will relate to the UV coupling of the theory. The integrable 
system has a spectral curve £ defined by the algebraic relation between x and another 
variable t: 

F(x, t) = det (T(x) - tl) = . (3.25) 

This curve is easily seen to have genus N — 1. The relation between the spectral curve 
and the integrable system is well known. The moduli space of the curve parameterizes 
the space of conserved quantities of the system. There are 3N variables subject to iV 
constraints 1)3.21)1 : hence, the dimension of the phase space is 2N. In particular, the 
action variables are given by integrating the 1-form xd log t around a set of basis 1-cycles 
on the curve. The conjugate angle variables are associated to a point in the Jacobian of 
S. 5 

The curve S has the form 

t 2 - 2tP(x) - h(h + 2)Q{x) = , (3.26) 

where 

2N 

Q(x) = Y[(x-m I ) (3.27) 



i=i 



with 

m 2 i-i = -Aj - ki , m 2i = -Xi + hi , (3.28) 

i — I, . . . , N, and where 

P( x ) = X N + H lX N - 1 + H 2 x N ~ 2 + ■ ■ ■ + H N . (3.29) 

Notice that the form of the curve is identical to the Seiberg-Witten curve of the M = 2 
theory with SU(N) gauge group and Nf = 2N fundamental hypermultiplets, and where 
h is expressed in terms of the coupling by [2U ESj 

%) = ^ , P.30) 



5 Note that the curve actually has genus N — 1 and so there is an apparent mismatch between the 
number of action-angle variables, N each, and the number of basis 1-cycles, or dimension of the Jacobian, 
of the curve, N — 1. We will resolve this conflict below. 
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where A(r) involves the automorphic function 




(3.31) 



and q = e mT . In particular, at weak coupling h = 32q + • • • . The parameters mj in ()3.27|) 
are not directly the bare masses, in fact 

2N 

mi = m b j + m b j/(2N) . (3.32) 

7=1 

In addition, the coupling r in the curve is not equal to the bare coupling of the theory 
[TH] . In general they are related by an instanton expansion of the form (jl.2j) . 

In (I3~2H|) the H^S®) basis for the Hamiltonians of the integrable system. Ex- 

plicitly, the first two Hamiltonians are 

i 

and 

H 2 = J2 (A^+^l+JOtsf Sp*-^ • (3-34) 

i<j 

3.3 Breaking to M = 1 

Now that we have established the precise relation between the twisted spin chain and the 
M = 2 field theory, we now consider the problem of breaking to Af = 1 using the integrable 
system. The procedure is well established [HI EI]. First of all, consider compactifying the 
theory to three dimensions on a circle of finite radius R. In the 3-dimensional effective 
theory, the gauge field in the [/(l)^" 1 unbroken gauge group can be "dualized" to scalar 
fields. These fields, along with the Wilson lines of the [/(l)^ -1 around the circle amass 
into complex scalar field. The N — 1 complex scalar fields are then naturally valued on 
a complex torus; in fact, precisely the Jacobian torus of S. In summary the Coulomb 
branch of the theory in three dimensions is identified with the complexified phase space 
of the integrable system and the split between the action-angle variables describes the 
moduli of the Coulomb branch of the four-dimensional theory plus the dual photons and 
Wilson lines. 



13 



Breaking to Af = 1, as in fjl.ip . is described by an effective superpotential on the 
Coulomb branch of the 3-dimensional theory. The resulting superpotential is simply, as 
described in the last section in detail, one of the Hamiltonians of the integrable system. 
From ()3.19j) . we have the following expression for the effective superpotential 



W int = Re Soo+ W(x)^ = Re Soo+ W(x) (- + ^j^O 
In particular, since W(x) contains the term — £x we have the constraint 



(3.35) 



{P> h{h + 2)Q' \ 
ReS - +X U + 2y(y + P) ) = °- (3 ' 36) 

This can be written 

h 2N 

Hi(S (j) ) = 2 = + X ^ h ■ ( 3 ' 3T ) 

i=i 

At the level of the integrable system experience with the finite U(N) k quiver theories [TO] , 
where the relative U(l) factor are also IR free, suggests that the correct way to freeze 
out the U(l) factor is via a Hamiltonian reduction. So as well as imposing ()3.37|) on the 
phase space, one also takes a quotient by the conjugate variable to Hi. Notice that H 1 
generates a simultaneous rotation of the spins around the S3 axes and so the con- 
jugate variable is this angle. In practice it is easier not to perform the quotient explicitly 
and as a consequence we shall find a degeneracy in our critical points corresponding to 
this unwanted degree-of-freedom. 



3.4 Example: SU(2) 

In this section, we consider the vacuum structure for the case of gauge group SU(2) and 
the quadratic mass deformation W(x) = fix 2 / 2 — £x. After imposing the constraint, and 
using ()3.35|) . we find an effective superpotential 

W ^ = 2(hTT) ' (3 ' 38) 

So vacua are given by critical point of Ei2- First of all, we solve the constraint ()3.37j) 
explicitly by choosing 

S® = y, S®=\ 1 + \ 2 -y. (3.39) 

and then impose the constraints (j3.21j) with Lagrange multipliers £j, i — 1,2. In order 
the that the four equations that result from the derivatives of H2 with respect to S^l 
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(2) 

and S{2, do not imply that the four aforementioned variables vanish (which would not 
be compatible with the constraints (|3.2ip ). we need 

66 + ^ = 0. (3.40) 

From this we find 

= - S [V +i (l + h)Si 2) ^ ^(i) = -»(1 + fr)£f } - 5f _ (3 41) 

The equation dH^jdy = can then be solved for y. This leaves the two constraints (j3.21|) 
which imply that the Lagrange multiplier £i satisfies a sixth-order polynomial equation 

- tekl£l + 128^ + 16(2(/i 2 + 2h - 2)k\ + 4A 2 + h(-2k\ + QXj + 2X 2 2 ) + 
h\-k\ + 2X 1 + 2X\M + 32h(h + 2)(-A; 2 + k 2 2 (l + />)(A 2 - A 2 ))£ 3 

+ 4/i(/i + 2)(-/i(/i + 2)k\ + 2((h 2 + 2h- 2)k\ + (h + 1)(2A 2 + h(X\ + A 2 ))))£ 2 

- 8h 2 {h + 2) 2 A; 2 ^ - h 3 (h + 2fk 2 2 = . 

Note that each of the 6 solutions is degenerate due to simultaneous rotations in the 
(Si, and (S[ 2 \ S^) planes. However, this is the expected and unphysical degeneracy 
in the conjugate variable to Hi. 

So the final result is that there are six vacua, which is the number expected (j2.14|) . 
For generic values of the parameters it is not possible to write down explicitly the critical 
value of the superpotential. However, let us consider their weak coupling, h — > limit. 
There are three pairs of solutions: 

(i) ft = 1 ± Xi/ki + 0(h) and f 2 = 0(h), which gives 



S {2) = (a,r } ,\ 1 + \ 2 ±k 1 }+0(h) , 



(3.43) 



where a 2 + if = k\ — (X\ + A 2 ± k 



\2 



(ii) ^ = 0(h) and f 2 = 1 ± X 2 /k 2 + 0(h), which gives 
S^=(t,r ] ,X 1 + X 2 ±k 2 ) )+0(h), 



(3.44) 
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where £ 2 + rf = k\ - (Ai + A 2 ± k 2 f . 



(iii) 6 = ±((kl - \l)h/2{\\ - fc2))V2 + o(/ l ) > w hich gives 

V 2 ^ 2 ^ ' (3.45) 
5( 2 ) = (a, 77, A 2 ) + 0(h) , 



where a 2 + if = A 



In all these solutions the degeneracy due to rotations in (£, 77) corresponds to the angle 
variable conjugate to Hi and is not physically relevant. Clearly the four vacua (i) and 
(ii) are the Higgs vacua with an expansion in instantons h ~ e mT at weak coupling. For 
these vacua we can write the first two terms in the weak coupling expansion as 



WU = ^,- ^'-^ + mj) h + 0^) (3.46) 



The two vacua (iii) are identified with the confining vacua since they have an expansion 
in terms of h 1 ^ 2 , in other words, finite instantons, at weak coupling. In these vacua 

W* nt = ±2/2(m 1 m 2 m 3 m 4 h/2) 1/2 + O(h) . (3.47) 

These results may be checked very directly by factorizing the Seiberg-Witten curve. 
Defining the curve as 

G(x, t)=t- 2P(x) - h(h + 2)Q(x)/t . (3.48) 

For the SU (2) case this genus one curve degenerates to genus zero when 

dG dG . A . 

These equations can be reduced to a single sixth-order polynomial equation for x: 

h(h + 2)Q'(x) 2 + AQ{x)P\x) 2 = , t = -h(h + 2)^^y ■ (3.50) 

The six solutions correspond to the four Higgs and two confining vacua and one can verify 
the result above. 
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4 The Holomorphic Matrix Model 



The effective superpotential of our supersymmetric field theory can also be calculated 
from a holomorphic matrix model in the way describing originally by Dijkgraaf and Vafa 
[21)] and then extended to include fundamental matter in [27] (see also the review j3j and 
references therein). 

The matrix model is defined by the partition function 

Z = J [d$] [dQ] [dQ]e-^ WtIc ^' m , (4.1) 

where W tre e(^, Q, Q) has the same form as the field theoretical tree level superpotential, 
but each $, Q and Q are now N x N, N x Nf and Nf x N matrices, respectively, not fields. 
Here N is a number of "colors" to be taken the large N with fixed physical quantities of 
N, Nf and S = g s N. The integral is to be understood in a holomorphic sense. 

The matrices Q and Q appear quadratically and can be integrated out. Gauge trans- 
formation can then be used to diagonalize $ at the expense of introducing a gauge- 
fixing determinant — the Vandermonde determinant. Denoting the eigenvalues of $ as Aj, 
% = 1, 2, . . . , N, the partition function becomes 

n d\ n^-(A, - Aj) h Ef=i w{Xi) 



We now take the limit N — > oo and g s — > with S = g s N fixed. In this limit, it is 
appropriate to introduce a density of eigenvalues 



PW = jJ26(\-K), (4.3) 

which is normalized as 

d\p(X) = l. (4.4) 
The free-energy of the model has the expansion 

-\xyZ = - 2 T x=2 + -T x=1 , (4.5) 
9s 9s 
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where 

T X=2 = S J d\p{\)W{\)- S 2 J d\d\' p{\)p{\')\og{\-\') , (4.6a) 

N f 

T x =x = SJ2 f d\ p(A) log(A - mj), (4.6b) 
1=1 

are a contribution from the planer diagrams (sphere) and diagram with one boundary, 
respectively. 

To leading order, the density of eigenvalues is determined from the saddle-point equa- 
tion 

2 / dX'f^j; = ±W(A) . (4.7) 

In the classical limit (S — > 0), the eigenvalues sit at one of the critical points of W(x). 
In other words, a classical configuration of eigenvalues can be described by the filling 
fraction Ni/N of the 2 th critical point. We are assuming that W(x) is a polynomial of 
degree n + 1. However, the saddle-point equation includes interaction terms coming from 
the Vandermonde determinant as well as from integrating out the matter matrices. At 
this point, we apply the usual hypothesis which states that the effect of these terms is 
to pull the eigenvalues away from the critical points along a set of open contours Cj, 
i — 1 . . . , n, in the neighbourhood of each critical point. The saddle-point equations (|4.7jl 
can be solved by introducing a resolvent 

uj(x) = [ dX-^- , (4.8) 
Jc x - A 

where C = U" =1 Cj. Note that u(x) is an analytic function over the complex plane with 
cuts along Cj. The density can be recovered from uj(x) as the discontinuity across the cuts 
C: 

uj{x + e) — uj{x — e) = —2iTip(x) , x G C , (4.9) 

where e is a suitable infinitesimal, so that x ± e lie just above and below the cut, respec- 
tively. In terms of u(x), the saddle-point equation becomes 

u(x + t) +u(x-e) = -W'(x) , xeC. (4.10) 

Since u>(x) goes to at infinity, the solution to the saddle-point equation ()4.10|) is imme- 
diate: 

2Su(x) = W\x) - ^W'(x) 2 + f{x) , (4.11) 
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where f(x) is an arbitrary polynomial of degree n — 1. This latter requirement ensures 

lim 2Su(x) = 0(l/x) . (4.12) 

It is useful to define 

y 2 = (W'(x) - 2Su(x)) 2 = W'{xf + f(x) . (4.13) 

The two sheets of y describe a hyper-elliptic Riemann surface S mm with genus — at least 
generically — of n — 1. 

The k constants in f(x) = J2]=i bjX^~ l are the moduli of the saddle-point solution, 
i.e. of the auxiliary surface E mm . These k moduli represent the freedom to choose the k 
filling fractions at each cut and we may change variables {bi} — > {Si}, where 

Si = g s Ni = I p{x)dx = —— "t / y(x)dx , (4.14) 



/ 1 ' 47TZ / a . 

where Aj is a contour which encloses the cut C,. In particular, 

S = y^Si = -y^— & y(x)dx= — (p y(x)dx = --^ , (4.15) 
where we pulled the contour off the cuts to the point at infinity oo + on the top sheet 



(y=y/W'(x) 2 + f(x)). 

We now have all the ingredients to define the glueball superpotential which is the 
effective superpotential of the physical field theory, in a particular classical vacuum, where 
the Si are interpreted as fields: 

Wg b = fa^ 1 + 7rzr5 <) + ■ (4 - 16) 

Here, r is the bare coupling of the field theory and the classical vacuum is specified by 
the Ni, Y^i=i Ni = N, the number of eigenvalues of the adjoint field $ lying at the i th 
critical point of W(x). 

To evaluate the derivative of T x =i in ()4.16|) . one can relate it to the variation of the 
free energy in bringing an eigenvalue in from infinity to the point at the end of the i th 
cut. Since the force on eigenvalue is y we have 
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The contribution from the matter fields, the third term in f)4.16|) . is easily evaluated, again 
one can show [2E| 

Nf / r°°~ \ 

T x=1 = | ( - / V dx + W(mj)) . (4.18) 
1=1 J m i 

In ()4.17|) and f)4.18|) . the integrals are actually divergent at infinity and need to be regu- 
larized. However, since we will work ultimately in the finite theory where the divergences 
cancel we shall not explicitly describe the details of the regularization process. In ()4.18jl . 
the point with x = mi can be either on the top or bottom sheet of £ and this freedom 
to choose is necessary in order to describe all the vacua of the physical theory [2HJ - Since 
we will later choose all the rrij on the bottom sheet when we engineer the Seiberg-Witten 
curve, we take the integral out to the point at infinity on the bottom sheet. If we need 
to specify whether x = mj is on the top or bottom sheet we write mf, respectively. 

Using ()4.17|) and ()4.18|) . the glueball superpotential ()4.16|) takes the form 

W gh = -J2^ y dx -jJ2f y dx -\Y, / ydx-W(mj) . (4.19) 

i=l J ^ i=l jA i 1=1 V m I J 

We now specialize to the case Nf = 2N and in this case the potential divergences at the 
infinite points cancel and the regular can be removed. 



4.1 Engineering the Seiberg-Witten curve 

The Seiberg-Witten curve itself can be extracted from the matrix model by choosing a 
suitably generic potential W(x) with order N + 1: 

JV 

W'(x)= f i N+1 H(x-( i ) . (4.20) 

i=i 

This potential acts as a probe for the M = 1 vacuum which can zero-in on any point on 
the Coulomb branch of the M = 2 theory. 

The matrix model curve S mm for an iV-cut solution at the critical point of the glueball 
superpotential, with one physical eigenvalue at each critical point Ni — 1, is then identified 
with the Seiberg-Witten curve where the {Q} are related to the moduli {a«} of the Af = 2 
Coulomb branch as in ()3.26|) . In addition the points x = mi are taken on the bottom 
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sheet, which we indicate as m 7 . We now proceed to find the critical point of the glueball 
superpotential in this case. 

It is important for what follows that the N — 1 quantities 

— ydx = dx , (4.21) 

dh y 2y K J 

for i — 1, . . . , N — 1, form a basis for the abelian differentials of the first kind on E. We 

will choose another basis for the abelian differentials {u>i}, i = 1 . . . , N — 1, normalized 

by 

/ ujj = k , i,j = l,...,N-l. (4.22) 

J Ai 

In addition, 

9 , 1 , 

—ydx = - Too+oo- , (4.23) 

dbi 2/i7v +1 

where on the right-and side we have t PiP2 , the meromorphic 1-form with simple poles at 
Pi and pi with residues =Fl, respectively. 

The glueball superpotential (|4.19J) (with Ni — 1 and Nf = 2N) can be written 

W ^ = "2 E M rfx - 47T-^ - 2 S / - ^ X " W(m/) ' (424) 

i=l , '- B i PV+l /=1 y-'m J y 

where Bi is the contour that goes from oo~ through the i th cut to oo + . To reach this 
expression we used the fact that J e °° ydx = | y dx. As we have already explained this 
expression must be regularized, although at the end in the Nf = 2N case the regulator 
can be removed. The careful treatment of the regularization process is given in 



Taking the derivative of (|4.24J) with respect to bi, i — 1, . . . , N — 2, (we will consider 
the b^-i equation below) and taking suitable linear combinations to change basis between 
{x l dx/y} and {uji}, we have 

N 2N oo- 2N 

= - J2 L - J2 / u * = ~ N / w * - Z) / Ui ' ( 4 - 25 ) 
j=i i=i i=i * m i 

where in the last line we have equality up to an element of the Jacobian lattice, pi + r^q^, 
Pi-, qi G Z. 

According to Abel's theorem, ()4.25|) implies that there exists a meromorphic function 
on E with a divisor (oo + ) (oo~) m^" • • -m^. This function is explicitly 

t = P(x)+y, (4.26) 
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Figure 2: The integral paths for the effective superpotential calculation. The curve has 
two sheets connecting N cuts and there are two distinguish infinite points, which we 
denote oo + and oo~. The integral cycles are defined in the figure. 

where P(x) = YliLi( x ~ a d * s a polynomial of degree N. In order that t has zeros at 
x = mi on the bottom sheet implies that 



P(x)+ 1 Q(x) = W'(xY + f(x) . 



(4.27) 



where Q(x) = Yli=i( x ~ m i)- These are 2N equations that determine {a^foj}. There 
is one remaining unknown 7 which will be fixed below. Hence, at the critical point the 
curve can be written 

y 2 = P(a;) 2 + 7 Q( x ) , (4.28) 

We recognize this as precisely the Seiberg-Witten curve (|3.26J) and in particular the {cij} 
are the moduli on the Coulomb branch. 



The remaining equation is obtained by taking the derivative with respect to b^. After 
using a Riemann bilinear relation with dlogt and d(y dx) / db^^i one obtains 



t(oo~ 



T 



_l lo l-y/T+7 

t(oo+) ni ° g 1 + ^/TTi ' 



(4.29) 



Note that this expression, for the case Nf < 2N, requires regularization, however, in the 
finite theory this is not necessary and a finite answer is obtained. Expressing 7 = h(h + 2), 
as in (|3.26|) . we find 

h = - - . (4.30) 



1 + e 71 
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Comparing with ()3.30|) . we see that the matrix model coupling, which we now denote as 
r mm , cannot be equal to the Seiberg-Witten coupling r sw directly, rather they are related 
by an expression of the form (jl.2j) . 

We can now see that the matrix model approach to the vacuum problem agrees pre- 
cisely with the integrable system treatment. Firstly the matrix model curve and Seiberg- 
Witten curve are identical, (j3~2fijl and IjOty . Secondly one can compare the critical value 
of the superpotential. Using once again a Riemann bilinear identity for dlogt and ydx 
one finds that the critical value of the glueball superpotential is 



a result that agrees precisely with W* nt in (j3.19|) . 

4.2 Effective superpotential for the quadratic potential 

In this section, we analyze the effective superpotential for the quadratic potential W(x) = 
fix 2 /2. This corresponds to looking at the one cut solution of the matrix model. Actually, 
in order to compare with the results of Section 13.21 we have to take into account the 
behaviour of the overall U(l) factor. The matrix model calculation on the face of it 
involves the gauge group U(N) rather than SU(N). However, the U(l) factor is, as we 
have described earlier, actually frozen out and is non-dynamical. However, in the matrix 
model approach the trace of <3> is non-vanishing. One finds 



In order to compare with our earlier results we can force the trace part of $ to vanish 
by considering, as previously, a more general potential W(x) = fix 2 /2 — £x. The extra 
parameter £ acts as a Lagrange multiplier which will be fixed by varying the effective 
superpotential. 

The second question that confronts us is how the one cut solution can describe the 
multiplicity of vacua that we described in Section EJ The origin for the multiplicity of 
vacua is solved by realizing that there is an ambiguity as to which sheet the points wi; lie. 
It is useful at this point to compare to brane configurations in Section 12.21 If mj lies on 
the bottom sheet then we interpret this to mean in brane language that the corresponding 




(4.31) 



r 




(4.32) 
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semi-infinite D4-brane ends on the right-hand NS5-brane, while, conversely, mj lies on the 
top sheet then the semi-infinite D4-brane ends on the left-hand NS5-brane. This means 
that the confining vacua, for example, have all the mi on the bottom sheet. 



For the one-cut solution, y 2 = (fix — £) 2 + 61 and the glueball superpotential is 

w gh = - 2 E ( / y dx - w ( m i)) , ( 4 - 33 ) 

where e\ is a point at the end of the cut. We can evaluate the integrals in this expression 
explicitly to give 

_ 2N 

W gh = ^"fri + ^E [fi 2 m 2 I -2£mi+fimif±(mi) + b 1 log (/im 7 +/ ± (m;)) -\b x log(-&i)J 

(4.34) 

where we have defined rhi = mi — £///. Note that f±(x) = ±a/ fi 2 x 2 + b\ depending upon 
whether x is on the top or bottom sheet, respectively. 

The critical point equation is 

2N 

II + /±(™/)) = (-bife mT . (4.35) 
7=1 

the critical value of the glueball superpotential is then 

^ 2N 

Kh = 7-Yl (V 2 ™ 2 i - 2£m 7 + /xmj/tCmj)) . (4.36) 
4/i 7=1 



4.3 The vacuum structure 

In Section we labelled the vacua by a subset {m^, . . . , mi r } C {mi, . . . , m 2 jv}. We will 
identify this vacuum with the matrix model saddle-point solution where {m^, . . . , mi r } 
are on the top sheet and the remainder mj, I <f_ {ii, . . . , I r }, are on the bottom sheet. In 
order to see why this is the correct identification, consider the classical limit r — > zoo. In 
this limit 61 is small and we can find a solution to (|4.35|) order-by-order in b\. To leading 
order we find 

t-i2N ~ 

(-61)^ = {2^ N -^e mr \\l= im l +■■■ . (4.37) 

U j= imj. 



24 



and 

The parameter £ is determined from dW* h /d£ = 0, giving to leading order 

For example, for the confining vacuum r = and all the m; are on the bottom sheet 
and to leading order £ = 0. In this case, to leading order 

W* gh = -fiN ( - e™ H m/J + • • • , (4.40) 

7=1 

which is the characteristic expansion in terms of fractional instantons. Correspondingly, 
in a Higgs vacuum, r — N — 1, so N — 1 of the mj are on the top sheet. In this case, 



N-l N-l TI^N ( _ , v^-'V-l 



* = ?( E m l + ( E ) + — ' N -i \ 2 e +■■■ ■ ( 4 - 41 ) 

Using the relation between the couplings in ()4.30|) . one can confirm that (|4.41|) and 
(|4.40|) for N = 2 match precisely with the results from the integrable system, ()3.46|) and 
(|3.47|) . respectively. 
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